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Abstract 

We report on the high-precision astrometric observations of maser sources around the Galactic Center in 
the SiO J = 1-0 v = l and 2 hues with the VLBA during 2001 - 2004. With phase-referencing interferometry 
referred to the radio continuum source Sgr A*, accurate positions of masers were obtained for three detected 
objects: IRS 10 EE (7 epochs), IRS 15NE (2 epochs), and SiO 6 (only 1 epoch). Because circumstellar 
masers of these objects were resolved into several components, proper motions for the maser sources were 
derived with several different methods. Combining our VLBA results with those of the previous VLA 
observations, we obtained the IRS lOEE proper motion of 76±3 km s~^ (at 8 kpc) to the south relative to 
Sgr A*. Almost null proper motion of this star in the east- west direction results in a net transverse motion 
of the infrared reference frame of about 30 ±9 km to the west relative to Sgr A*. The proper-motion 
data also suggests that IRS lOEE is an astrometric binary with an unseen massive companion. 
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1. Introduction 

It is believed that the center of the Milky Way harbors 
a super massive black hole (Lynden-Bell & Rees 1971). 
Indeed, interferometric observations confirmed the pres- 
ence of a compact radio source there (Balick & Brown 
1974; Ekers et al. 1975; Lo et al. 1975), and the unre- 
solved source was named Sgr A* (Brown 1982). In the last 
decade, dramatic improvements in resolution of infrared 
instruments enabled to detect stellar proper motions un- 
der gravitational influence of the central massive object, 
which indicates an enclosed mass of ~ 3.7 x 10^ Mq (Ghez 
et al. 2005; Eisenhauer et al. 2005). An upper limit of 1.8 
km s^^ (2(7) to the Sgr A* random motion in the extra- 
galactic reference frame also suggests that it must contain 
at least a mass of 0.4 x 10*^ Mq (Reid & Brunthaler 2004). 
These findings provide compelling evidence that Sgr A* 
is a source associated with the super massive black hole. 
Because of large extinction at optical wavelengths, obser- 



vations have been attempted at X-ray, infrared, and radio 
wavelengths (Muno et al. 2005; Rafclski et al. 2007; Shcn 
et al. 2005). 

Collocation of the positions of the Galactic Center at 
different wavelengths is a key issue for the research. The 
position of the compact radio source (Sgr A*) appeared 
to be coincident with an extended 2.2 //m source observed 
within ± 1" (Becklin & Neugebauer 1975). To improve the 
accuracy of the position of Sgr A*, Menten et al. (1997) 
carried out accurate radio measurements of the positions 
of SiO and H2O masers in red giants relative to the con- 
tinuum emission of Sgr A*; the alignment of the radio 
and A'-band (2.2 ^m) infrared reference frames achieved 
an accuracy of 0.03". More precise radio/infrared reg- 
istration is important for determining relative positions 
between Sgr A* and the IR flare (Nayakshin 2005; Eckart 
et al 2006), or between Sgr A* and periapses of stellar 
orbits (Mouawad et al. 2005); in fact, the periapsis of SO- 
16 was estimated to be about 45 AU (6 milliarcsec [mas]) 
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away from Sgr A* (Ghez et al. 2005). Therefore, attempts 
were made toward this direction with radio and infrared 
telescopes (Reid et al. 2003; Reid et al. 2007), reaching an 
accuracy of the concordance down to nearly one mas. 

An obstacle for position measurements in radio is a 
temporal variation of spatial distribution of circumstcl- 
lar masers. The brightness distribution of maser emission 
is not uniform in the circumstellar shell, but rather it is 
random and patchy. SiO masers in a typical mira are 
emitted at a few stellar radii, i.e., at about 2 x 10^^ cm 
(Cotton et al. 2004). This corresponds to 1.6 mas at the 
distance of the Galactic Center (8 kpc). Therefore a po- 
sition of the central star estimated from a single maser 
feature may have considerable uncertainty. Due to its 
high sensitivity, past measurements of proper motions of 
SiO masers around Sgr A* were made mainly with the 
VLA (Very Large Array); a synthesized beam of 80 mas 
X 40 mas towards Sgr A* at 43 GHz, though they advo- 
cated "combined" VLA and VLBA (Very Long Baseline 
Array) results for only a single-epoch VLBA observation 
(Reid et al. 2007). Though the VLBA is a very pow- 
erful tool to investigate a spatial distribution of maser 
features (a synthesized beam of 1.8 mas x 0.5 mas at 43 
GHz), it is limited by the sensitivity (Sjouwerman et al. 
1998). In the case of the SiO masers in the Galactic cen- 
ter sources, peak line intensities of the SiO J = 1-0 v = 1 
and 2 transitions are about 0.3 Jy or less. This is a se- 
vere restriction for the VLBI observations, which require 
high S /N ratios of the target maser lines in normal data 
calibration procedures. However, Deguchi et al. (2002) 
found that SiO maser emission from IRS lOEE occasion- 
ally flares; the intensity reached 1.5 Jy in the J = 1-0 
V = 1 line in April-May 2000, a factor of 5 stronger than 
the normal intensity (sec Izumiura et al. 1998). Because 
infrared monitoring observations demonstrated a periodic 
light variation of this object (Wood et al. 1998), it is rea- 
sonably expected that the SiO maser brightening occurs in 
every few year. Therefore, the positional accuracy of IRS 
10 EE could greatly be improved if the VLBI observations 
are made at the time of SiO flares. 

In this work, we demonstrate abilities and limitations of 
the VLBA to measure proper motions of weak SiO maser 
sources near Sgr A*. We reports on our few- years attempt 
to improve the positional accuracies of IRS 10 EE and two 
other SiO maser sources in the Galactic center. The rela- 
tive positions of IRS lOEE SiO masers with respect to Sgr 
A* were obtained with unprecedented accuracy. However, 
because maser features have been spatially resolved with 
the VLBA, refined analyses are necessary to extract the 
position of the central star from the data. Combining 
them with the infrared proper motions of this star, we de- 
rive the relative velocity of Sgr A* to the infrared reference 
frame, which gives a constraint on dynamical motions of 
various objects surrounding the Galactic Center. 

2. Observations and data reduction 

The VLBI observations of SiO masers around Sgr A* 
were conducted between 2001 May and 2004 March (8 
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Fig. 1. Visibility amplitude plotted versus M-j; distance. Tlie 
solid curve shows the best fitted visibility curves for a circular 
Gaussian with a zero baseline flux of 0.927 Jy and FWHM of 
0.712 mas, which is consistent with the size measured at 43 
GHz by Lo et al. (1998). 

epochs in total), using the Very Long Baseline Array 
(VLBA), a facility of the National Radio Astronomy 
Observatory (NRAO). The present dataset consists of 
three series of observations, which are summarized in 
Table 1. Our primary target stars were chosen from table 
2 of Reid et al. (2003). They are located within 50" from 
Sgr A* (R.A. = 17:45:40.0392, Decl. = -29:00:28.151 in 
J2000). In this paper, we followed the naming of SiO 
sources adopted by Reid et al. (2003). 

The first series of observations was made in 2001. The 
system temperatures (including atmospheric noise) were 
100-300 K. Received signals were sampled at a rate of 8 
M sample per second per base band channel (BBC) with 
2-bit quantization in only left-hand circular polarization, 
and recorded at the rate of 128 Mbps (16M bit x 8 
BBCs) using the VLBA terminal. In data correlation, five 
different processings were performed, in which the track- 
ing centers were set at Sgr A*, IRS lOEE, IRS 7, SiO 6, 
and SiO 11. To observe both v = 1 and v = 2 lines, we di- 
vided eight BBCs into two groups; one group at 42.820587 
GHz and another at 43.122027 GHz to cover the v = 2, 
J = 1-0 and u = 1, J = 1-0 transitions of SiO. Each BBC 
band was divided into 64 spectral channels. Each band 
yielded a spectral resolution of 62.5 kHz, equivalent to 
0.435 km s^^ in velocity, and covered a velocity range of 
27.84 km s^^. The four 4 MHz bands of one group cen- 
tered on velocities of —120, —60, —20, and 60 km s~^ with 
respect to the local standard of rest (LSR). In this first 
series, a correlator accumulation time was set to 1.05 s 
to cover a field of view of 6.2" from the tracking centers. 
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Table 1. Summary of VLBA observations. 



Session nick- Date v BW Aiy Acc. period detected object 

name name (yy/mm/dd) (MHz) (kHz) (s) 



A 146a 2001/05/16 1,2 4x8 62.5 1.05 IRS lOEE 

B 146b 2001/07/31 1,2 4x8 62.5 1.05 IRS lOEE 

C 146c 2001/10/05 1,2 4x8 62.5 1.05 failed 

D 177a 2003/01/11 1,2 8x4 125 0.13 IRS lOEE 

E 177b 2003/03/12 1,2 8x4 125 0.13 IRS lOEE, IRS 15NE, SiO 6 

F 177c 2003/10/06 1,2 8x4 125 0.13 IRS lOEE 

G 131b 2004/03/08 1 16x4 62.5 0.39 IRS lOEE, IRS 15NE 

H 131c 2004/03/20 1 16x4 62.5 0.39 IRS lOEE, IRS 15NE 



Note. — V is the vibrational state of the SiO J — 1-0 transition involved in the observation. Acc. period means an 
accumulation period. 



Therefore, not only the star located at the tracking cen- 
ter but also any other stars were observable within a field 
of view. For example, IRS 15 NE was within the field of 
view in which IRS 7 was used as a tracking center. The 
separation between IRS 15NE and IRS 7 (nearest to IRS 
15NE among five stars used as tracking center) is 5.9". 
Also, the separations between IRS 28 and IRS lOEE, and 
between IRS 12N and Sgr A* are 10.4" and 7.6", respec- 
tively. Therefore, although IRS 28 and IRS 12N were with 
in the same band, significant losses of sensitivity were ex- 
pected. 

The second series of observation was made in 2003. The 
weather conditions were good on 2003 Jan 11 and Mar 3, 
but the PT, KP, LA, OV and FD stations were in bad 
weather on October 6. To expand the field of view, the 
correlator accumulation time was set to be 0.13 s, yielding 
the field of view became about 50" from Sgr A*. Though 
the main purpose of this observation was to detect stellar 
masers as many as possible, the number of available BBCs 
were restricted. We used four BBCs of 8MHz each into 
two groups; one for v = 1 and another for v = 2. The spec- 
tral band was divided into 64 channels, yielding a spectral 
resolution of 125 kHz, equivalent to 0.86 km s~^ in veloc- 
ity, and covered a velocity range of 55.04 km s~^. The 
two BBCs in one group were centered at LSR velocities of 
—60 and —10 km s~^. As a result, we were able to observe 
IRS lOEE, IRS 15NE, SiO 6, SiO 11, IRS 17, and SiO 12. 

The third series of observations was made on 2004 Mar 
8 and March 20 to investigate the structural variation of 
Sgr A*. After the observations, three different correlation 
processings were performed; the tracking centers were set 
at Sgr A*, IRS lOEE, and IRS 9. Two BBCs of a band- 
width of 16 MHz each were used, both of which were set 
to the v = l line. Each BBC was divided into 256 spectral 
channels, giving a spectral resolution of 62.5 kHz, equiv- 
alent to 0.43 km s~^ in velocity, and it covered a velocity 
range of 110.08 km s^^. The two BBC bands were cen- 
tered at LSR velocities of —6 and 339 km s~^. In this 
series, the correlator accumulation time was set to 0.39 
s, yielding a field of view of 16.7". Therefore, not only 



stars at the tracking center but also IRS 15NE and IRS 
28 were included within several fields of view. The sep- 
arations between IRS 15NE and IRS lOEE and between 
IRS 28 and Sgr A* were 8.8" and 10.9 ", respectively. As 
a result, we were able to observe IRS lOEE, IRS 9, IRS 
15NE and IRS 28. Table 1 gives a detail of the setup for 
these observations. 

At all of the epochs, the continuum source NRAO 530 
was observed approximately every one hour with an inte- 
gration time of seven minutes to check the system perfor- 
mance and to calibrate the clock parameters. 

The data reduction was made in the standard manner 
with the Astronomical Image Processing System (AIPS). 
Initial delay and phase calibrations were made using the 
calibrator NRAO 530. Single-band delays were deter- 
mined with fringe-fitting and then applied to the Sgr A* 
visibilities. Accurate delays and rates were determined 
with global fringe-fitting directly to the Sgr A* visibili- 
ties. Sgr A* was detected on the baselines at u-v distances 
smaller than about 200 MA. Figure 1 shows the visibil- 
ity amplitude versus u-v distance for Sgr A* observed on 
2003 March 12. The data was averaged in each single 
frequency band (one channel) after residual delay, rate, 
phase and amplitude corrections. Bad data due to spuri- 
ous amplitude and unstable phase were flagged out using 
the AIPS task UVFLG. Then, the hybrid imaging pro- 
cess, which involved iterative imaging with CLEAN and 
self-calibration, was made. We obtained antenna-based 
solutions for correcting the atmospheric phase fluctuation. 
These calibration solutions were applied to the all maser 
objects. This method is known as the "phase-referencing 
technique", and makes it possible to determine the rel- 
ative positions of maser features with respect to Sgr A*. 
Differential aberrations were compulsorily corrected in the 
procedure UVFIX. After the phase referencing calibra- 
tion, Doppler correction of target maser data was applied 
using the AIPS task CVEL. Then we created image cubes 
of the z; = 1 and v = 2 masers using the AIPS task IMAGR. 
The initial images contained 512 x 512 pixels with a pixel 
size of 0.2 mas x 0.2 mas, and a taper was applied to the 
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Table 2. Position offsets of detected SiO masers from Sgr A*. 
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Jy km B-1) 


(mJy km B-l) 




( km s~^) 
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IRS lOEE 












A 


1 


1 


7684.43 


0.08 


4209.21 


0.25 


252 


47 


2 


-27.9 


0.24 


0.20 


A 


2 


1 


7684.40 


0.04 


4209.33 


0.10 


556 


63 


4 


-28.3 


0.08 


0.24 


B 


1 


1 


7684.55 


0.05 


4208.85 


0.17 


199 


39 


4 


-28.7 


0.05 


0.09 


B 


2 


1 


7684.32 


0.04 


4208.68 


0.14 


108 


16 


9 


-28.8 


0.48 


0.52 


D 


1 


1 


7683.55 


0.06 


4204.30 


0.12 


373 


73 


4 


-27.9 


0.30 


0.40 


D 


1 


2 


7681.22 


0.06 


4208.67 


0.19 


470 


80 


4 


-27.0 


1.31 


0.56 


D 


1 


3 


7680.88 


0.09 


4209.07 


0.17 


673 


103 


3 


-27.9 


0.22 


0.35 


D 


2 


1 


7683.24 


0.09 


4204.18 


0.14 


773 


107 


7 


-27.9 


0.24 


0.27 


D 


2 


2 


7681.20 


0.09 


4209.22 


0.24 


511 


99 


3 


-27.9 


0.16 


0.07 


D 


2 


3 


7680.70 


0.10 


4209.96 


0.22 


547 


103 


1 


-27.0 


0.20 


0.02 


E 


1 


1 


7682.08 


0.02 


4205.77 


0.07 


416 


41 


4 


-27.9 


0.05 


0.09 


E 


2 


1 


7683.66 


0.05 


4205.00 


0.14 


314 


52 


5 


-28.8 


0.18 


0.22 


E 


2 


2 


7682.20 


0.03 


4205.71 


0.09 


525 


53 


4 


-27.9 


0.05 


0.08 


F 


2 


1 


7683.81 


0.04 


4205.68 


0.13 


261 


39 


4 


-27.0 


0.34 


1.31 


G 


1 


1 


7684.26 


0.02 


4203.17 


0.05 


1231 


87 


13 


-28.7 


0.12 


0.12 


G 


1 


2 


7682.82 


0.04 


4203.84 


0.12 


570 


84 


7 


-27.9 


0.09 


0.29 


H 


1 


1 


7684.34 


0.03 


4203.29 


0.06 


138 


16 


6 


-28.7 


0.06 


0.06 














IRS 15NE 












E 


1 


1 


1217.24 


0.03 


11291.46 


0.09 


506 


60 


8 


-15.5 


0.70 


0.32 


E 


2 


1 


1217.36 


0.05 


11291.38 


0.09 


774 


83 


9 


-17.3 


0.35 


0.57 


G 


1 


1 


1216.35 


0.06 


11286.67 


0.15 


565 


96 


5 


-14.2 


0.10 


0.23 


H 


1 


1 


1216.37 


0.06 


11286.19 


0.15 
SiO 6 


370 


69 


3 


-14.2 


0.04 


0.54 


E 


1 


1 


35234.92 


0.03 


30662.85 


0.09 


542 


62 


5 


52.0 


0.04 


0.09 


E 


2 


1 


35234.83 


0.06 


30662.37 


0.15 


389 


69 


4 


52.0 


0.92 


0.37 



Note — V is the vibrational state of SiO J — 1-0 involved in the observation. 
Nc is the number of components detected in channel maps. 
Visr is the velocity of the strongest component. 

da and ds are standard deviations of the R.A. and Decl. offsets, respectively, 
of the components in channel maps. 



data, which reduced the spatial resolution and made the 
initial search for emission more efficient. The final image 
contained 512 x 512 pixels with a pixel size of 0.05 mas 
X 0.05 mas. 

The spatial and velocity structures of maser emission 
were more or less varying from epoch to epoch. The peak 
positions and flux densities of all maser components were 
derived from the elliptical Gaussian fitting to the compo- 
nents using the AIPS task JMFIT. If more than one com- 
ponents exist at same velocity channel, we estimated the 
position of each component by a series of single-Gaussian 
component fits. Also, if the structure is too complicated 
in a channel map, we consider an area of a peak flux as a 
single Gaussian component, and estimated the position. 

3. Results and discussions 

3.1. Distributions of maser features around stars in the 
Galactic Center 

SiO masers in the Galactic Center have low flux densi- 
ties, typically 60-500 mJy. Moreover, they are spatially 
resolved out in the long baselines of the VLBA. Therefore, 
it results in a low detection rate in the present observa- 



tions. Only three sources, IRS lOEE, IRS 15NE, and SiO 
6 (a new detection with the VLBA), were detected with 
the shortest 10 baselines at the signal-to- noise ratio (S/N) 
above 4. Figure 2 shows integrated intensity maps of IRS 
lOEE at three epochs and Figure 3 shows those of IRS 
15NE and SiO 6 in the SiO J = 1-0 w = 1 and 2 lines. 
The integrated-intensity map covers a velocity range of 

- 8 km s-i centered at Msr 27.0 km s"^ for IRS lOEE, 

a range of '--^ 7 km s~^ centered at V\sr ~ —13.3 km s^^ for 
IRS 15 NE, and a range of ~ 8 km s~^ centered at V\si ~ 
52.0 km s~^ for SiO 6. These VLBA maps show several 
maser features distributed in areas with a diameter of 3 
mas with velocities ranging a few km s~^. 

In IRS lOEE (=OH359.946— 0.048), the SiO maser 
emission was extended in an area of 2 mas x4 mas (Figure 
2), though apparent elongation to the north-south is an 
artifact by the elongated beam. It was resolved into sev- 
eral features. We simultaneously observed two transitions 
(SiO J = 1-0 v = l and 2) except the last two epochs. The 
location of the v = I emission well coincided with that of 
u = 2. Furthermore, strong time variation was found. The 
SiO masers flared on 2003 January (session D). The to- 
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Fig. 2. Integrated intensity maps of the SiO J = 1-0 = 1 and 2 lines for IRS lOEE in July 2001 (top), March 2003 (middle), and 
March 2004 (bottom). The v = 2 line was not observed at the last epoch. The half-power beam is shown at the lower left of each 
panel. Contours are drawn at every noise-level step above 3 times the noise levels, 20.7 (top left), 25.4 (top right), 18.6 (middle left), 
21.9 (middle right), and 26.3 mjy km s~^ beam~^ (bottom left), except for the middle-left and bottom panels, which use a step of 
twice the noise-level. 
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Fig. 3. Same as figure 2, but for IRS 15 (upper panels) and SiO 6 (lower panels) in Jan 2003. Contour levels are drawn at every 
noise-level step above 3 times the noise levels, 24.7 (top left), 31.0 (top right), 22.3 (bottom left), and 23.2 mjy km s~^ beam~^ 
(bottom right), except for the bottom left panel which uses a step of twice the noise levels. 



tal intensity reached about 1.6 Jy km s~^, a factor of few 
stronger than the normal intensity level. The masers were 
split into several features in an extended area of about 5 
mas X 3 mas; the north-west features were not present 
at the other epochs (see Table 2). Because of the large 
spread of emission at this epoch, we discarded the data 
of this session from the proper motion fitting described in 
the next section. 

A scale of 1 mas corresponds to 1.2 x 10^"* cm at a dis- 
tance of 8 kpc. For comparison, a typical mira variable, 
TX Cam, exhibits a ring structure of SiO maser emission 
with a diameter of ~ 2 x 10"'^'* cm at a distance of 450 
pc (28 mas; Diamond et al. 1994), which is comparable 
with the extent of the SiO masers in IRS lOEE. In the 
case of a supcrgiant, for example, for VY CMa, the dis- 
tribution of SiO maser features is extended in an area of 
80 mas (~ 120 AU) at the distance of 1.5 kpc from the 



Sun (Miyoshi et al. 2003; Shibata et al. 2004), so that it 
would be extended in an area of about 10 mas if placed 
at the Galactic center. The spectral type of IRS lOEE 
was classified to MOIII (Figer et al. 2003) and photomet- 
ric observations suggested that IRS lOEE is a long period 
variable (Tamura et al. 1996; Wood et al. 1998; Blum et 
al. 2003; Peeples et al. 2007). The small diameter of SiO 
maser emission of this star is consistent with these previ- 
ous findings. 

For the case of the IRS 15NE masers, the SiO emis- 
sion was slightly resolved and extended in an area with a 
diameter of about 2 mas (upper two panels of Figure 3). 
However, only a single maser feature appeared both at two 
epochs. This object was also identified to the AGB star 
(Genzel et al. 1996; Genzel et al. 2000), which is consistent 
with a measured small maser emission size. However, the 
IR spectrum of this star is blended by radiation from a 
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massive hot star (Najarro et al. 1997, or see a note in ta- 
ble 2 of Blum et al. 2003), and has been classified to WN8 
(Paumard et al. 2006; Martins et al. 2007). Therefore, an 
attention must be paid for the astrometric use of this star. 

The SiO 6 masers (lower two panels of Figure 3) were 
also partially resolved over the area of about 2 mas. The 
extent of the SiO masers suggests that the object is an 
AGB star. OH and H2O masers have been found in 
this star (Levine et al. 1995; Yusef-Zadeh & Mehringer 
1995; Sjouwerman & van Langevelde 1996; Lindqvist et 
al. 1992a). Though Levine et al. (1995) and Yusef-Zadeh 
& Mehringer (1995) considered that the H2O emission 
was associated with a supergiant or a molecular cloud 
near Sgr A*, subsequent observations (Sjouwerman & van 
Langevelde 1996) revealed that the H2O maser is asso- 
ciated with the OH/IR star (OH359.956— 0.050) with a 
typical doubly peaked spectrum with a velocity separa- 
tion of 22 km s~^. The light variation period of this star 
is 607 d ([GMC2001] 3-5 or V4928 Sgr; Glass et al. 
2001; Peeples et al. 2007). The SiO masers have previously 
been detected several times at Nobeyama (Izumiura ct al. 
1998; Miyazaki et al. 2001; Deguchi et al. 2002; Deguchi 
et al. 2004), and intensities of SiO masers seem to not 
vary drastically. Our result also confirmed that this is not 
a supergiant. 

3.2. Proper motions of the maser stars 

Stellar masers at the Galactic center are quite faint. 
Therefore, we used the in-beam phase-reference technique 
together with a long integration time. Noise levels in in- 
dividual channel maps reached ^ 20 mJy per beam. In 
order to measure accurate proper motions of individual 
SiO maser stars, we tried to track the central star posi- 
tion of masers with several different methods. The first 
one is to take the average position of all the maser fea- 
tures detected with S/N > 5 in the integrated intensity 
map in each epoch. These components arc listed in Table 
2. Though formal uncertainties of positions for individual 
components are quite small (^ 0.1 mas), the data points 
are scattered in a scale of a few mas. Therefore, it is nec- 
essary to re-evaluate the magnitudes of uncertainties for 
obtaining the position of the central star. If the maser 
features appear randomly in all directions referred to the 
central star, the average position must be close to the 
center position of the star. We estimated a standard de- 
viation of the position offset of a maser component from 
the central star to be about 0.8 mas (here we assumed the 
SiO maser feature appears randomly in a flat disk with a 
radius of 1.6 mas (3 stellar radii for a typical AGB star). 
Therefore, we added this value to the formal position un- 
certainty of the individual component. Figure 4 shows the 
linear fits to the R.A. and Deck offsets from Sgr A* as 
a function of time using both J ~ 1-0 v = 1 and 2 lines 
of IRS lOEE. The results obtained with this method are 
shown in the first row of table 3. 

The second method is to track a single maser feature. 
For this method to work, we have to select a suitably 
strong maser feature, which must last during the obser- 
vational period. For this purpose, we applied the criteria 




2 3 
YR - 2000 (yr ) 




01 2 3 4 S 

YR - 2000 (yr) 

Fig. 4. Linear fit to IRS lOEE positions relative to Sgr A* in 
riglit ascension (upper panel) and in declination (lower panel). 
Filled and unfilled circles denote the SiO J = 1—0 v = 1 and 2 
lines, respectively. The error bar (1.5 stellar radii; 0.8 mas) is 
also indicated in each data point. The data of session D were 
discarded. 

such as signal-to-noise ratio higher than five in the in- 
dividual channel map, detections in the several adjacent 
channel maps at the same position, and detections at three 
or more epochs at the same velocity. With these selection 
criteria, we determined the center position of the —27.9 
km s~^ component in sessions B, C, E, and G in both 
of the V ~ 1 and 2 lines for the IRS lOEE masers, and 
computed a proper motion by fitting a straight line to the 
positions. The results are listed in the 2nd row of Table 
3. 

The uncertainties in the second row of Table 3 arc quite 
small compared with those in the first row because we 
used the formal measurement uncertainties of individual- 
component positions in the second method. Because a 
single maser feature is expected to move relative to the 
central star by an amount of about 20 km s~^, we need to 
correct the motion by about 0.52 mas yr~^ at most for the 
central star. Therefore, taking into account this effect, we 
can state that the proper motions obtained by the second 
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method are consistent with the values obtained by the first 
method, provided that the —27.9 km s^^ component is an 
clump at the south-west edge of the expanding envelope 
of IRS lOEE, though this does not clearly reveal from the 
data. 

From OH scattering-size measurements at 1.6 GHz (van 
Langevelde et al. 1992; Frail et al. 1994), angular broaden- 
ing of a maser feature by interstellar scattering at 43 GHz 
is expected to be up to about 1 mas at the Galactic cen- 
ter. Therefore, the center of a single Gaussian component 
is an average position of several maser features within the 
same velocity channel, making a situation worse, even if 
we track the same single velocity component. Bower et 
al. (2001) suggested that the center position of Sgr A* 
shifts with frequency by 2-5 mas to the east between 8.4 
and 5 GHz, which was attributed to a difference of source 
structure between two frequencies (in Sgr A* and position 
reference sources). However, this can also be attributed 
to the frequency-dependent refraction of radio waves in 
plasma toward the Galactic center. This refraction ef- 
fect could cause position offsets between radio and optical 
measurements by about 0.1 mas at maximum at 43 GHz, 
which is not negligibly small. However, it does not influ- 
ence the proper motion measurement unless it systemati- 
cally varies during the observational time period. 

The third row of Table 3 gives the IRS lOEE proper 
motions obtained from Table 2 of Reid et al. (2007). 
They reasonably agree with the present results. Note that 
the uncertainties involved in the first method are slightly 
larger than those of Reid et al. (2007). This is because 
we used the position uncertainty of 0.8 mas for deducing 
the central star position. Figure 5 shows a comparison 
of our measured positions (averages) and those of Reid et 
al. (2007). Note that the proper motion of IRS lOEE in 
the east-west direction is negative in our measurements, 
while that of Reid et al. (2007) is slighty positive; this can 
be clearly seen in Figure 5. 

If we combine both results, we expect to obtain better 
values of proper motions for these stars. For our data to 
be compatible with those of Reid et al. (2007), we took all 
the maser components with high S/N (> 5) listed in Table 
2 and adopt an uncertainty of 0.8 mas. which was used in 
the first method. The combined results are shown at the 
4th, 8th and 10th rows in Table 3. Note that the com- 
bined results of our VLB A and Reid et al.'s data (mostly 
with the VLA), do not necessarily give the middle value, 
but somewhat different one, because the proper motion 
is a time derivative between positions. The position un- 
certainties, which are used to weight on the data points, 
influence the derived values; here we also added the same 
uncertainty 0.8 mas mentioned above to the measurement 
uncertainties of Reid et al. (2007). The obtained posi- 
tions of IRS lOEE relative to Sgr A* in unit of arcsec are 
written 

Aa = 7.68390[±0.00029]" - 0.000057[±0.000097]" 

X (y [year] -2000.0), (1) 
AS = 4.21198[±0.00029]" - 0.002007[±0.000093]" 
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Fig. 5. IRS lOEE positions relative to Sgr A*. Filled cir- 
cle and unfilled triangle indicate the average positions of IRS 
lOEE SiO masers obtained at each epoch in this work (with 
VLBA), and those observed by Reid et al. (2007) with VLA 
(except 1996.41 with VLBA), respectively. The number in- 
dicates the epoch of individual observations. The error bar 
of the data is also indicated. The dotted line shows the best 
linear fit for combined results. 

X (r[j/ear] -2000.0). (2) 

The combined result offers currently best values of the 
proper motions of IRS lOEE and we use these values in 
the later discussion. Additionally, we noticed in figure 5 
that there still remains a small periodic shift of the star 
position from the best fit line; to the east during 2000- 
2001, and to the west during 2003, returning back to the 
cast in 2004-2006. This shift seems to be a motion due 
to a binary system. We discuss on this issue in Appendix 
1. The analysis demands further corrections in proper 
motions, but they are negligibly small. Therefore, we dis- 
carded them in the later discussions. For IRS 15NE, we 
get 

Aa = 1.22364[±0.00039]"-0.001873[±0.000104]" 

X (r[year] -2000.0), (3) 

A(5 = 11.31044[±0.00041]"-0.005761[±0.000111]" 

X (y[year] -2000.0), (4) 

We can compare these results with nominal NIR posi- 
tions of IRS lOEE and IRS 15NE, and pin-point the Sgr 
A* positions on their NIR images. However, because the 
difference between Reid et al. (2007) and the combined 
result is quite small (^0.1 mas yr^^) compared with the 
position uncertainties of stars in the infrared (typically 
more than 10 mas), we cannot effectively improve the Sgr 
A* position in the high resolution NIR images better than 
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Table 3. Proper Motions of IRS lOEE, IRS 15NE, and SiO 6. 
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H"y 


<7y 


note 








(km s^i 


) fnias vr~^ 


) fmas vr""^") fi 




^) fnias vr "'^l 




IRS lOEE 


6 


1,2 


all 


-0.32 


0.21 


-1.96 


0.21 


this work {SIN > 5) 




4 


1,2 


-27.9 


-0.43 


0.02 


-2.25 


0.07 


this work (a feature) 




6 


1 


-27 


0.04 


0.07 


-2.09 


0.07 


Reid et al. (2007) 




12 


1,2 


all 


-0.06 


0.09 


-2.01 


0.09 


combined 


IRS 15NE 


3 


1,2 


all 


-0.92 


0.79 


-4.90 


0.79 


this work [SjN > 5) 




3 


1 


-13.7 


-1.60 


0.05 


-5.70 


0.06 


this work (a feature) 




6 


1 


-12 


-1.96 


0.07 


-5.68 


0.12 


Reid et al. (2007) 




8 


1,2 


all 


-1.87 


0.10 


-5.76 


0.11 


combined 


SiO 6 


4 


1 


52 


2.58 


0.43 


1.99 


0.52 


Reid et al. (2007) 




5 


1,2 


all 


2.58 


0.26 


1.70 


0.44 


combined 



Note. — Ne is the number of epochs involved. 



Reid et al. (2007) did. 

The proper motions of several hundred stars have been 
measured on the NIR images. They are individual star 
motions relative to the average motion of the ~ 400 "well- 
behaved" stars within ~ 30" of the central cusp (Reid et al. 
2003; Reid et al. 2007). Referring to the IRS lOEE proper 
motions, we can derive the motion of Sgr A* relative to 
this NIR reference frame, i.e., relative to the central cusp. 
From the IR proper motions (/ij,^ and /i^^ in table 4 of 
Reid et al. 2007), we can compute the relative motion as 
(A^^, Afiy)= (0.79 ±0.24, 0.08 ±0.27) [mas yr^i] using 
IRS lOEE, and (-0.53 ±0.48, -0.53 ± 0.37) [mas yr"!] 
using IRS 15NE. The proper motion in R.A. direction 
obtained from IRSIOEE is significantly larger than the 
standard errors (3.3 a significance). The same motions 
obtained using IRS15NE involve larger errors, so that they 
are discarded. We conclude that the relative motion of 
the Sgr A* to the nuclear star cluster is 0.79 mas yr^^ 
to the east (or 29.9 ±9.1 km s~^ at a distance of 8 kpc). 
Reid et al. (2007) also derived the RA relative motion of 
24 ±9 km s~^ from their data (2.7 a), which is consistent 
with this work. It should be noted, however, that the 
value strongly depends on the accuracy of proper motion 
measurements of stars on the NIR images and sampling 
of "well behaved" stars. It is supposed that blending and 
binarity of reference stars in the dense central star cluster 
may strongly affect to the results, especially the motion 
of IRS lOEE in this work. Therefore, above conclusion 
is somewhat provisional, and future NIR high-precision 
astromctry can examine the validity of above result. 

3.3. A speculation on the anomalous motion of Sgr A* 

The average of the velocities of ^ 400 stars in the cen- 
tral cusp was estimated in an accuracy of about 5 km s~^ 
(a statistical deviation due to random motions of the ref- 
erence stars; Reid et al. 2007), which is much larger than 
a brownian motion of the central black hole in the central 
cusp (^^0.2 km s ^^; Merritt et al. 2007). Therefore, above 
mentioned anomalous motion of Sgr A* with respect to 
the NIR reference frame of about 30 km is consider- 



ably larger than the statistical deviation, and needs to be 
explained. Here, we give a speculation on this anomalous 
motion. 

Suppose that a majority of stars in the Galactic cen- 
ter rotate around the barycenter of the Galaxy (here, 
Sgr A*) in the same direction. If a sampling is com- 
plete, the average transverse velocity of these stars with 
respect the barycenter would be zero due to point sym- 
metry. However, if thick dust is present near the Galactic 
center, magnitude-limited sampling of stars is biased to- 
wards the front-side stars of the Galactic center, because 
the stars behind are fainter. In fact, Davidge (1998) and 
Alexander (1999) studied K-band extinction within a few 
arcminutes of Sgr A* and found that the increase of ex- 
tinction by about AAk ~ 0.3 toward Sgr A* (down to 
nearly 30" radius). A simple calculation based on the 
cusp density law leads about 39% of stars seen towards 
Sgr A* within a projected distance of 1.5 pc 40") to 
be actually located at radii larger than 1.5 pc from Sgr 
A*. Therefore, the sampling of these distant stars are 
possibly biased to the near side due to extinction by thick 
dust located near Sgr A*. If we assume a circular rota- 
tion of stars around Sgr A* in the same sense, this bias 
results in the preferential selection of the stars moving in 
one direction. An estimation of the magnitude of this ef- 
fect gives a net motion of about 29 km s^^ [« 120 km s~^ 
xO.39 X 0.5/(0.61 + 0.5 x 0.39)] ; here all stars behind Sgr 
A* at radii over 1.5 pc arc assumed to be not seen and the 
average velocity of stars within 1 . 5 pc is zero) . Though the 
value is strongly dependent on the geometry and velocity 
distribution of the stream, it yields an acceptable magni- 
tude of net motion of NIR reference frame with respect 
to the barycenter of the Galaxy. The observed direction 
of the motion of the NIR reference frame is to the west 
referred to Sgr A*. Therefore, the direction of the stellar 
stream has an angle of about 60 degree from the Galactic 
plane and the sense of rotation is opposite to the Galactic 
rotation. 

It has been discussed that an intermediate-mass black 
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hole (IMBH; 10^ - lO^M©) may strongly influence on 
the dynamics of the nuclear star cluster containing a super 
massive black hole (Portegies Zwart et al. 2006; Berukofi^ 
& Hansen 2006; Matsubayashi et al. 2007). Numerical 
simulations found that the IMBH alters eccentricities of 
stellar orbits and, as a result, alters the power-law index 
of the central cusp considerably. IMBH may also be re- 
sponsible for ejecting hypervelocity stars from the central 
star cluster (Baumgardt et al. 2006). Eiscnhauer et al. 
(2005) and Paumard et al. (2006) found two young-star 
disks with total masses of ~ 1.5 x 10^ Mq and ages of 
about 6 Myr; one rotates clockwise, and another counter- 
clockwise. Therefore, dynamical structure of stars at the 
inner few arcscc seems considerably complex. Because the 
400-600 IR position reference stars, which were used for 
estimating the average motion of the nuclear star cluster, 
are mostly located out of Sgr A* by more than a few arc- 
sec, they are not much affected by the IMBH which is 
supposedly at a few arcsec from the Galactic center (e.g., 
IRS 13E; Schodcl et al. 2005). Moreover, the observed 
power-law index of the cusp density between 4" and 100" 
is approximately 2 (for example, see figure 7 of Genzel et 
al. 2003). This is close to the steady-state power low in- 
dex of 7/4 of the cusp star density with a super massive 
black hole (e.g., Bahcall & Wolf 1976), indicating that the 
IMBH docs not influence to the motion of stars outside a 
few arcsecond. Therefore, it is unlikely that the streaming 
motion associated with the NIR reference frame is caused 
by IMBH. The total mass within a radius of 15 pc is about 
3 X 10'' Mq (for example, see Figure 6 of Deguchi et al. 
2004), which is by about one order of magnitude larger 
than the mass of the central black hole. Therefore, the 
streaming motion of stars discussed here must be a dy- 
namical property related to the central core of the inner 
Galactic bulge. 

4. Summary 

The 8-epoch VLBA observations of maser sources in the 
Galactic center were made in the SiO J — 1-0 v = 1 and 2 
lines at -43 GHz during 2001 May - 2004 March. The SiO 
masers in IRS lOEE had been detected through all of the 
epochs, while those in IRS 15NE and SiO 6 were detected 
only at several epochs. The radio continuum source Sgr 
A* were also observed within the same single-dish tele- 
scope beam. Accurate astrometric positions of the SiO 
maser features relative to Sgr A* were obtained at each 
epoch with the simultaneous VLBI imaging of the maser - 
continuum source pairs. However, because of the spatial 
distribution of maser features in the circumstellar enve- 
lope, the derivation of the accurate values of proper mo- 
tions for the maser sources was not straightforward. The 
proper motions were derived with two different methods: 
one, using the average positions for all the detected fea- 
tures with high signal-to-noise ratios, and the other, using 
the positions of long- lasting features; both methods gave 
very similar values. Thus, the precise measurements of 
proper motions of maser sources relative to Sgr A* were 
made. Combining our VLBA positional data with the for- 
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Fig. 6. Orbit of IRS lOEE. The filled circles in the line graph 
indicate annually averaged positions of SiO masers after linear 
motions are subtracted. The best-fit ellipse is shown in solid 
curve. 

mer measurements by VLA, and with NIR proper motion 
of IRS lOEE, we found an unexplained net motion of Sgr 
A* relative to the nuclear star cluster of about 30 ± 9 km 
s-i. 
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Appendix 1. IRS lOEE binary system 

We noticed that the positions of IRSIOEE in Figure 
5 deviate systematically from the regression line to the 
east during 2000-2001, to the west during 2003, and re- 
turning back to the east in 2006. It looks like a wavy 
motion of an astrometric binary star with an unseen com- 
panion. Therefore, it is worth to pursue a possibility of 
the binary system with a period of about 6 yr. In fact, 
Peeples et al. (2007) found that the X-ray source CXOGC 
J174540. 7-290024, which was fisted in table 2 (on-line 
material) of Muno et al. (2003), is located within 0.8" of 
IRS lOEE (their IRS 10*), and concluded that IRS lOEE 
is a binary system involving a white dwarf. This source 
has X-ray luminosity Lx = 1.1 x 10"^^ erg s~^. which is 
considerably strong for a AGB star with a white dwarf 
companion 10^^ erg s"-'^). 

Therefore, it is interesting to attribute the residuals 
from the best-fit line [equations (1) and (2)] to a binary 
motion of IRS lOEE, and to check whether or not the de- 
rived orbital parameters of the binary system are reason- 
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Table 4. Orbital elements for the IRS lOEE binary system 



Element 


unit 


value 


Aao 


mas 


0.137 (±0.186) 


Mo 


mas 


-0.115 (±0.033) 




mas yr~^ 


-0.029 (±0.065) 


Afly 


mas yr~^ 


0.013 (±0.012) 


T 


yr 


5.65 (± 0.12) 




mas 


0.73 (±0.41) 


h 


mas 


0.35 (±0.10) 


ei 




0.87 (±0.24) 


i 


degree 


« (±4) 



Note. — Aao J and A6o are the corrections for the center 
positions of the ellipse, and A/x^ and A^y are corrections 
for the proper motions, relative to the quantities given 
in equations (1) and (2). Quantities, T, ai, bi, e, and i 
are the period, semi-major and semi-minor diameters of 
the primary, ellipticity, inclination of the orbit to the sky 
plane, respectively. 

able. However, a strong restriction for a binary model is 
obtained from the radial velocity of IRS lOEE; the radial 

velocity of SiO mascrs of IRS lOEE {Visr 27 km s'^ 

) did not change more than a few km s~^ in the past 10 
years. This fact gives a severe condition on the binary or- 
bit; it must be almost face-on to us. Note that this binary 
model is valid only when the apparent motion revealed 
in figure 5 is attributed to a binary motion. However, we 
cannot completely deny the other possibility that it is due 
to a systematic periodic variation in SiO maser emission 
distribution in the circumstellar envelope of IRS lOEE. 

Though the procedure obtaining orbital elements of a 
binary system has been well established (for example, 
see a brief introduction by Asada et al. 2007), it is not 
straight forward in the present case because of large po- 
sitional errors of the central star involved in SiO maser 
astrometry. The best-fit line obtained in section 3.2 does 
not represent the trajectory of the barycenter of the bi- 
nary system, but merely an approximate average position 
of the primary (here IRS lOEE) of the binary system. It 
has to be corrected for the elliptical orbit of the primary. 
In order to obtain the corrections for the motion of the 
center of the ellipse of the primary, we again fit the R.A. 
and Dec. residuals of IRS lOEE from the average motion 
[equations (1) and (2)] by introducing the linear shift plus 
a sin curve (as a function of time), and determined the 
period, corrections for the ellipse center, and corrections 
for proper motions. After correcting these quantities for 
the center of the ellipse, we plotted the residual Xc and 
Yc in Figure 6, and fit an ellipse to these points using the 
formula of Asada et al. (2007). Then, we obtained the 
orbital elements for the primary, which are listed in Table 
4. Because of the above mentioned restriction on the or- 
bit being face-on (i.e., i w 0), the ellipse we obtained on 
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Fig. 7. Radial velocity variation of IRS lOEE SiO masers. 
The filled and unfilled circles indicate the observed peak ve- 
locities of the SiO J = 1—0 v = I and 2 lines, respectively, in 
the present paper and in the past literatures. The thick curve 
indicates the velocity variation expected from the elliptical 
orbit given in Table 4 with an inclination angle of j = 0.7°. 
The broken curve with smaller amplitude is a sine-curve fit 
for comparison. 

the sky plane must be very close to the real binary orbit. 
Alternatively, the inclination angle i can be determined by 
the condition of constant areal velocity (the Kepler's sec- 
ond law; Asada et al. 2007). However, such a procedure 
did not give a reasonable inclination angle in the present 
case because of large scatter of the data points. Figure 7 
shows the observed radial velocities of SiO emission peak 
of IRS lOEE in the literatures (e.g., see Izumiura et al. 
1998; Deguchi et al. 2002) and in the present paper. 
The very small variation of radial velocities (AVisr ^ 3 km 
s""'^) seems to be fit by the elliptic orbit model given in 
table 4 with an inclination angle of i = 0.7°, though such 
a small variation of peak velocities can also be attributed 
to a variation of velocity fields in the expanding envelope 
without a binary model. 

From the obtained orbital elements, we can compute 
an approximate mass of the unseen companion using the 
Kepler's third law. 



mi + J712 = /T^ 



(5) 



where mi and m2 are masses of primary and a compan- 
ion (in unit of solar mass) and T is the orbital period (in 
unit of year), and a (in unit of astronomical unit) is the 
mean separation of the binary, which can be written with 
the semi-major diameters of elliptic orbits of two stars as 
a = fli +02. Here, they are constrained by ai mi = 02 m2. 
In order to determine the mass of the companion, we have 
to give a mass of the primary (in addition to ai and T) . 
The mass of IRS lOEE (primary here) is deduced to be 3 

- 6 Mq, because it is a relatively massive AGB star. From 
equation (5), we obtain the mass of the companion to be 
in a range of 13.1 - 16.2 Mq. Because this massive com- 
panion is not seen on infrared images, it cannot be a main 
sequence star; it must be a black hole (a white dwarf or 
neutron star is excluded because of the large mass). The 
mean separation between two stars is in the range of 1.1 

- 1.2 X 10^^ cm, and the mean orbital speed of IRS lOEE 
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is about 30 km s^^. The periastron distance between two 
stars is in the range of 1.8 - 1.9 x 10^^ cm, indicating 
that the secondary grazes the extended atmosphere of the 
primary (re ~ 1.5 x 10^'^ cm for MOIII) at the periastron 
passage, when a drastic event may happen. However, the 
periastron distance strongly depends on the uncertainty 
of 1 — e of the orbit, and it may vary by a factor of 3. 

Estimated values of the parameters for the IRS lOEE bi- 
nary system seem to stay within tolerable ranges, though 
they are slightly extreme in a sample of stellar-mass black 
hole binaries (McClintock & Remillard 2003). Therefore, 
under the condition that the observed nearly null radial 
velocity variation of SiO masers indicates the face-on ge- 
ometry of the orbit, we conclude that IRS lOEE is an 
astrometric binary possibly with a black hole companion. 
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